Induction kinetics of a conditional pH stress response system in Escherichia coli 
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The analysis of stress response systems in microorganisms can reveal molecular strategies for 
regulatory control and adaptation. Here, we focus on the Cad module, a subsystem of E. colics 
response to acidic stress, which is conditionally activated at low pH only when lysine is available. 
When expressed, the Cad system counteracts the elevated concentration by converting lysine 
to cadaverine under the consumption of H^, and exporting cadaverine in exchange for external 
lysine. Surprisingly, the cad operon displays a transient response, even when the conditions for 
its induction persist. To quantitatively characterize the regulation of the Cad module, we have 
experimentally recorded and theoretically modeled the dynamics of important system variables. 
We establish a quantitative model that adequately describes and predicts the transient expression 
behavior for various initial conditions. Our quantitative analysis of the Cad system supports a 
negative feedback by external cadaverine as the origin of the transient response. Furthermore, the 
analysis puts causal constraints on the precise mechanism of signal transduction via the regulatory 
protein CadC. 
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Introduction 

During their natural life cycle, gastrointestinal bacte- 
ria are faced with acid stress while passing the extreme 
low pH of the stomach and being exposed to volatile 
fatty acids in the intestine. Escherichia coWs remark- 
able ability to sustain growth over multiple decades of 
concentrations [T] and its potential to survive ex- 
tremely low pH is implemented by a battery of pH home- 
ostasis [2 -4 and acid tolerance systems [5,-^. In recent 
years, it was increasingly recognized that each of these 
subsystems is specifically activated under certain envi- 
ronmental conditions [6, 8 , while the orchestration of 
the different responses is just beginning to be explored. 
However, a system-level study of the acid stress response 
requires a detailed quantitative analysis of the individual 
modules. 

One of the conditional stress response modules is the 
Cad system [9UT3]. which is induced only when acidic 
stress occurs in a lysine-rich environment. The three 
principal components of the Cad system are the enzyme 
CadA, the transport protein CadB, and the regulatory 
protein CadC, see Fig. [l] The decarboxylase CadA con- 
verts the amino acid lysine into cadaverine, a reaction 
which effectively consumes H+ [14 . The antiporter CadB 
imports the substrate, lysine, and exports the product, 
cadaverine. Together, CadA and CadB reduce the in- 
tracellular H+ concentration and thereby contribute to 
pH homeostasis |15j. The cytoplasmic membrane 
protein CadC not only senses the external conditions 
[121 mi Hi] 5 but also regulates the response by binding 
directly to the DNA and activating the transcription of 



cadBA [18 . Like other members of the ToxR family [19], 
CadC thereby performs signal transduction in a single 
component, without the phosphorylation step employed 
by two-component systems [20 . Fig.[l]also depicts the ly- 
sine permease LysP, which is not part of the cad operon, 
but essential for its function, since CadC senses lysine 
indirectly via interaction with LysP [121 IISl [III EI] • In 
contrast, the external (periplasmic) pH is believed to be 
sensed directly by CadC, through a pH-dependent con- 
formational transition |T7l [22]. The signal integration 
performed by CadC then assures that CadA and CadB 
are produced only under the appropriate external condi- 
tions of low pH and lysine abundance. 

However, CadC also senses a third input, which seems 
surprising from a physiological point of view: External 
cadaverine binds to CadC [21] and represses the long- 
term expression of the cad operon [T2[. As cadaverine 
is the end product of the decarboxylase reaction, it was 
suggested that it accumulates in the medium and causes 
a delayed transcriptional down-regulation of cadBA ex- 
pression Although many stress response systems 
display a similar transient response [23l [24], their reg- 
ulation strategy appears to be fundamentally different: 
For instance, the osmo- stress response of yeast directly 
follows its stimulus (low osmolarity) and remains active 
until the osmolarity returns back to physiological levels. 

In this study, we explore whether a negative feedback 
via external cadaverine can account for the transient re- 
sponse of the Cad system on a quantitative level. In par- 
ticular, we are interested in the factors that determine 
the duration and the amplitude of the transient response, 
and ask how the addition of external cadaverine affects 
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FIG. 1: Qualitative Model of the Cad System in E. coli (sim- 
plified). The Cad system is conditionally activated by low 
pH and high lysine levels. Lysine inhibits the repressive ef- 
fects of the lysine permease LysP on the receptor CadC, while 
low pH activates CadC directly. The active form of CadC 
activates transcription of cadBA, encoding the lysine decar- 
boxylase Cad A and the lysine/cadaverine antiporter CadB. 
The Cad system imports lysine, decarboxylates it under con- 
sumption of a cytoplasmic proton and exports the product 
cadaverine in exchange for another lysine molecule. The net 
effect of these reactions is the expulsion of a proton from the 
cytoplasm. Finally, it is believed that external cadaverine 
deactivates CadC. In our experiments we recorded the time- 
evolution of the variables in grey boxes. 



these characteristic quantities. It is known for instance, 
that external cadaverine reduces the long-term activity of 
the Cad system [12 , but one would like to know whether 
it shortens the duration of the transient expression pulse 
or reduces the amplitude of the pulse. Or does it affect 
both of these properties? 

To address these questions, we quantitatively mea- 
sured the dynamics of the Cad system in three important 
variables at high time resolution: the cadBA transcript, 
the activity of the lysine decarboxylase CadA, and the 
concentration of excreted cadaverine. Based on the ex- 
isting qualitative model, cf. Fig. [T] we formulate a quan- 
titative model for the Cad system, and test its agree- 
ment with the experimental response dynamics. We find 
that our quantitative model coherently describes the dy- 
namical response of the wild-type Cad system within 
a physiological parameter regime. The available data 
constrains the key biochemical parameters to a narrow 
regime. For instance, we infer the effective in vivo de- 
activation threshold for the Cad system and compare 
it to a previously measured in vitro binding threshold 
[21 . Using the quantitative model, we formulate pre- 
dictions for the response dynamics of the Cad system 
under conditions with initially added cadaverine, and in 
a mutant strain with a defective lysine permease LysP. 
The successful experimental validation of these predic- 
tions strongly supports the existence of the postulated 
feedback inhibition mechanism via cadaverine in the Cad 



system. Finally, we discuss the causal constraints of our 
results on the signal transduction mechanism by CadC, 
helping to discriminate between two contradicting mod- 
els. 



Results 

Transient expression dynamics 

To probe the transient response of the Cad module, we 
first grew E. coli (strain MG1665) to exponential phase 
at pH 7.6 in minimal medium. We then induced the Cad 
module by transferring cells into fresh minimal medium 
with 10 mM lysine and buffered at pH 5.8. The induction 
defined the starting point, t = min, for our measure- 
ments of the response, which we performed initially at 
intervals of 5 minutes, until t = 30 min, and then at 
longer intervals of 30 minutes. To quantify the response, 
we assayed the cadBA mRNA level, the specific CadA 
activity, and the external cadaverine concentration, see 
Materials and Methods for all experimental details. The 
results are shown in Fig.|2] Over the entire 4 hour period 
of the experiment, the external pH, shown in Fig. |2](a), 
remained low, even slightly decreasing from the induc- 
tion level. Transcription of the cad operon began imme- 
diately after induction, and mRNA rapidly accumulated 
as shown in Fig. [2] (b). At t ^ 25 min the mRNA level 
peaked and then rapidly decreased, reaching its low pre- 
induction level at about t = 90 min. The response on 
the protein level, as quantified by the specific activity 
of CadA shown in Fig. [2] (c), was slower, exhibiting a 
slight delay after induction and reaching a plateau level 
at t ~ 60 min, which was sustained over the time of the 
experiment. The activity of the Cad module led to the 
production and secretion of cadaverine, which accumu- 
lated in the medium as shown in Fig. [2](d). 

The transient expression of the Cad module shown 
in Fig. |2] is in qualitative agreement with previous in- 
duction experiments which studied the system in a less 
quantitative manner [13]. The biochemical mechanism 
for the transient behavior remained unclear, however. It 
was suggested that the external cadaverine level exerts a 
negative feedback on the activity of the regulator CadC 
[m [13]. Alternatively, the Cad module might, for in- 
stance, directly affect and control the level of its input 
stimuli. In other words, the activity of the Cad module 
might reduce the external lysine concentration below its 
induction threshold, or shift the external pH level outside 
its range for induction. 

Dose-response curves 

To address these possible alternative explanations, and 
to characterize the ranges and the intensity of the total 
response under our experimental conditions, we next de- 
termined the "dose-response" behavior of the Cad mod- 
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FIG. 2: Induction kinetics of the cad operon in E. coli 
MG1655. (a) The Cad system was induced at t=Omin by 
a shift from pH 7.6 to pH 5.8 and simultaneous addition of 
10 mM lysine. The time evolution of the cadBA mRNA (b), 
the CadA activity (c) and the extracellular cadaverine con- 
centration (d) was determined as described in Materials and 
Methods. All values are average values from duplicate repeti- 
tions. For a detailed discussion see main text. 



ule. We have seen above, see Fig. |2](c), that the CadA 
activity reaches a steady-state plateau about 60 min after 
induction. We take this plateau value as a proxy for the 
total (cumulative) response of the Cad system and study 
its dependence on the input signals. To this end, we in- 
duced the Cad module with different external pH levels 
and initial lysine concentrations, and assayed samples at 
least 90 min after induction for their CadA activity, see 
Materials and Methods for details. Fig. [s] (a) and (b) 
show the lysine- and pH-dependence of the response, re- 
spectively. The data in Fig. [s] (a) indicates that when 
induced with pH 5.8, the Cad module is barely active at 
lysine concentrations below 0.5 mM, whereas it is fully 
active for lysine levels exceeding 5 mM. In between these 
values, the activity increases sigmoidally with the induc- 
ing lysine concentration. Similarly, at a given lysine in- 
duction level of 10 mM, the activity depends sigmoidally 
on the inducing pH, see Fig. [3] (b), with no significant 
activity above pH 6.8, and full activity at pH 5.8 and 
below. 

Taken together, the pH-dependence of the total re- 
sponse shown in Fig.[3](b) and the time series in Fig. [5] (a) 
shows that the pH level did not leave the range for in- 
duction during the course of our experiment in Fig. |2j 
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FIG. 3: Dose-response curves of the wild-type Cad system, 
(a) Dependence of the long-term CadA response on the in- 
ducing lysine concentration. The specific CadA activity was 
determined 8 h after induction with pH 5.8 and stated ly- 
sine concentration, see Materials and Methods. Similarly, the 
dependence on the inducing pH in (b) was obtained by in- 
duction with 10 mM lysine and stated pH. Here activity was 
determined 90 min after induction. The solid lines show the 
fit result of our quantitative model. 



Hence, the transient behavior of the cadBA expression 
is clearly not mediated by a decrease of the external pH 
stimulus. 

The dose-response curves of Fig. [3] characterize the 
input-output behavior of the system when the Cad mod- 
ule is regarded as a "black box" signal processing unit. 
In particular, we can read off the apparent activation 
thresholds of the Cad module, i.e., the pH and lysine 
levels at which the module displays half-maximal activ- 
ity. Such apparent thresholds constitute the first level of 
description in a top-down system analysis. Conversely, in 
a bottom-up analysis, the first level of description is via 
biochemical interaction parameters, while effective pa- 
rameters such as apparent thresholds emerge from the in- 
terplay of molecular interactions. Quantitative modeling 
and analysis of this interplay is the only way to connect 
between the two levels of description. In the following, 
we want to make such a connection, and then leverage it 
to estimate the molecular activation thresholds of CadC 
from the apparent behavior of the Cad module. 



Construction of a quantitative model 

It is clear from the above that a minimal quantitative 
model of the Cad module must describe the integration 
of the input signals pH and lysine, as well as the effect 
of cadaverine on the activity of CadC. Furthermore it 
must describe the regulation and expression of the cadBA 
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operon, and the functioning of the CadA and CadB pro- 
teins, such that we may relate the quantitative model 
to the observed dynamics of the three system variables 
monitored in our experiments. Our construction of such 
a quantitative model is guided by the qualitative model 
in Fig. [T] and the known biochemistry of the Cad module. 

Signal integration. The three external signals known 
to affect the activity of the Cad module are the time- 
dependent lysine concentration, /(t), cadaverine concen- 
tration, c(t), and the pH. The membrane protein CadC, 
which receives and combines these signals into a single 
response [12] [161 [13 ^ is constitutively expressed ? 
hence we take the total amount of CadC per cell, Co, to 
be constant. The signals then modulate only the fraction 
of active CadC molecules per cell, C{t)/Co. We assume 
that the different signals regulate CadC independently, 
such that the CadC activity is described by the product 
form 



C(t)/Co = /(pH(t)) X gm) X h{c{t)) 



(1) 



Indeed, all experimental data available show no indica- 
tion for a coupled effect of the input signals on CadC 
[ini [21] . Another assumption implicitly made by Eq. ([T]) 
is that the fraction of active receptors is always equili- 
brated to the current levels of the input signals, i.e., it 
does not depend on the signal levels prior to the time 
t. This assumption is also plausible, since the typical 
timescale for conformational transitions in receptors (see, 
e.g., Ref. [25]) is much shorter than the timescale of our 
experiments. The functions /, g and h in Eq. ([T]) take 
on values between and 1, and are assumed to be of the 
Hill form typical for cooperative binding reactions. The 
pH-dependence is parameterized as 



/(pH) = 



1 



1 + 10 



pHp 



(2) 



with pHo denoting the pH value at which / reaches half- 
maximal activity and ApH determining the width of the 
response curve. Similarly, the lysine and cadaverine de- 
pendence take the form 



5(0 
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(3) 



where Ki and Kc are the effective in vivo activation 
thresholds for the direct and indirect regulatory interac- 
tions of lysine and cadaverine with CadC. As usual, the 
Hill coefficients, ni and nc, parameterize the cooperativ- 
ity of the binding reactions and determine the maximal 
sensitivity for signal detection. The difference in form 
between g{l) and h{c) stems from the fact that lysine 
activates the Cad module whereas cadaverine represses. 
On the other hand, the difference to Eq. ([2| is because 
the pH is logarithmically related to the H+ concentration. 

Transcriptional regulation. In its activated conforma- 
tion, CadC directly binds to the cadBA promoter and ac- 
tivates cadBA expression [18 . Generally, transcriptional 



regulation in bacteria can be described by quantitative 
"thermodynamic" models [26] [27]. For the present case, 
an appropriate form for the transcriptional activity from 
the Pcad promoter as a function of the abundance of ac- 
tive CadC, is derived in Materials and Methods. 
The resulting rate equation for the time evolution of the 
mRNA level m then takes the form 



dt 



m{t) = Vrr 



l^{C(t)/Kc?f 
l^{C{t)/Kcf 



\m m{t) , (4) 



with the basal transcription rate the degradation rate 
A^, the fold-change / between basal and maximal tran- 
scription rate, and Kc denoting the binding threshold 
for CadC-DNA binding. The particular choice of the ex- 
ponents in the first term is motivated by the observation 
that the cadBA promoter appears to be regulated by two 
binding sites for dimeric forms of CadC [18 . 

Kinetics of enzyme expression and catalysis. On the 
protein level, we have a similar interplay of synthesis and 
decay as in Eq. 0, 

j^A{t) = iypm{t)-XpA{t), (5) 

where A is the abundance of CadA per cell, and Vp and 
Xp are the translation and degradation rate, respectively. 
The level of the transporter CadB is taken to be propor- 
tional to that of CadA. Since they are translated from 
the same mRNA, this means that we neglect possible 
post-transcriptional regulation, for which there seems to 
be no experimental indication [TT]. We also assume that 
we can subsume the transport and turnover of lysine to 
cadaverine through CadB and CadA by a single effective 
reaction, since little is known about the microscopic rates 
and affinities of the coupled transport and decarboxylase 
reactions. As detailed in Materials and Methods^ this 
assumption leads us to 



lit) 



(6) 



This simplified reaction corresponds to an effective 
Michaelis-Menten process with external lysine as the sub- 
strate, an effective maximal lysine turnover rate Vjjiax-, 
and an effective Michaelis constant i^^. Implicit in 
Eq. ([6| is also the assumption that growth of the bac- 
terial population over the time period of the experiment 
is negligible. For the external cadaverine level, we assume 
the flux balance 



d , . <^ w X 



(7) 



dt dt 

implying that the sum of external lysine and cadaver- 
ine is conserved at all times, l(t) + c(t) = const. This 
flux balance appears justified, given experimental results 
with a LysP-deficient mutant strain which we report and 
discuss further below. We take the pH to be a constant 
over the duration of our kinetic experiments, since the 
pH changes only very little in our buffered medium, and 
the Cad module is not very sensitive to the pH over this 
regime (see above). 
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Data interpretation with the quantitative model 

We will now demonstrate that the simple quantitative 
model constructed above is indeed a powerful tool. First, 
we test to what extent this model is compatible with the 
data sets reported above. For this test, we also include 
another data set from Neely et al. [12], who determined 
the cadaverine-dependent dose-response of the Cad mod- 
ule, see Table |Tj This additional data further constrains 
our model and probes the consistency with the existing 
literature. 

In total, the quantitative model has 14 parameters. We 
constrained each of these to a range inferred from typical 
physiological values and other information in the litera- 
ture, see Table [TTj We then fitted our model to all data 
sets simultaneously using standard least-squares mini- 
mization of the residual ^is described in Materials 
and Methods. The curves corresponding to the best fit 
parameters are shown in Figs, [s] and [4] (blue lines). The 
overall agreement with the experimental data is good, 
both for the response dynamics in Fig. |4] (a)-(c) as well 
as for the dose-response in Fig. [3] and Table [ll Note how- 
ever, that not all model parameters are individually well 
constrained by the data. This becomes apparent by plot- 
ting the correlations between the quality of the fit, char- 
acterized by the residual and the fit parameters, as 
shown in Fig. [s] As the fit becomes better (lower x^)^ 
most parameter values are confined to a narrow inter- 
val, indicating that the information contained in the ex- 
perimental data accurately determines their values. For 
instance, the cadBA mRNA lifetime determines the 
decay time of the transient expression peak in Fig. |4](a), 
and is therefore strongly constrained in our model. In 
contrast, some of the parameters display a wide varia- 
tion even at the lowest x^ values, e.g. the transcription 
and translation rates, v^n and Vp. For these cases, where 
individual parameters are 'sloppy' [28], certain combina- 
tions of these parameter are well constrained by the data 
sets. Pairwise scatterplots, as shown in Fig. SI of the 
Supplementary Material^ identify correlations and anti- 



TABLE I: Repression of the long-term Cad response by 
cadaverine (data taken from Neely et al. [E])- The 
cadBA expression {central column)^ as determined from the 
/3-galactosidase activity of a cadA — lacZ fusion, was mea- 
sured in cells that were grown in medium pH 5.8 with 10 mM 
lysine and the indicated cadaverine concentrations for 3 h [12] . 
The model values in the right column show the fit result of 
our quantitative model. 
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FIG. 4: Fit of our quantitative model (blue lines, left column) 
to the experimental data of the wild-type induction kinetics 
{blue circles). The parameter- free prediction of our model for 
a LysP-deficient mutant (red lines, right column) displays rea- 
sonable agreement with the experimental data (red triangles). 
The dashed and solid black lines in (e) and (j) show the lysine 
and cadaverine contribution {g{l{t)) and h{c(t)) in Eq. ^) to 
the signal integration function in Eq. ([Tj). 



correlations between the parameters and help to reveal 
the appropriate combinations. For instance, the product 
of the transcription and translation rates is much better 
determined by the data than the individual rates. 

It is noteworthy that the best-fit value for the effective 
Hill coefficient Uc for the regulation of CadC by cadav- 
erine is close to 3, and is relatively well constrained by 
the data. This suggests that a molecular mechanism for 
cooperativity is at work, possibly a multimerization of 
CadC proteins in the membrane. Another interesting 
observation from Table Ull concerns the half-life of cadBA 
mRNA, which was well-constrained by the data to a value 
of almost 14 min. A global analysis of RNA half-lifes in 
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TABLE II: Parameters of the quantitative model. 

LB - lower bound; UB - upper bound; Estimated parameter values are shown as (best fit value) =b ^_ , where cr^ and cf~ indicate 
the asymmetric standard errors in positive and negative direction, respectively, see Eq. (14) in Materials and Methods. 
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pHo 
ApH 



Sensory module 

Threshold for CadC activation by lysine 

Threshold for CadC inactivation by cadaverine 

Hill coefficient for CadC regulation by lysine 

Hill coefficient for CadC regulation by cadaverine 

pH threshold for CadC activation 

width of the transition from active to inactive CadC 



Expression module 

Co/K(- Total CadC per cell i 
binding 



relation to the threshold for CadC-promoter 



basal transcription rate 

fold-change between basal and maximal transcription rate 
effective translation rate 
mRNA half-life (= In 2/ Am) 
protein half-life (= In2/Ap) 

maximal rate for lysine turnover via CadA and CadB 

effective Michaelis constant for lysine turnover via CadA and CadB 
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estimated from data in Fig.p|(b) 
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100 


2,±% 


mM 
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Escherichia coli [29 found an extremely short half-life of 
less than 2 min for the cadBA mRNA, suggesting that an 
active degradation mechanism is involved. Our Northern 
blot data for the lysP211 mutant, shown in Fig. [4] (f), 
does indeed suggest a rapid decay of the mRNA at high 
levels directly after the peak, followed by a slower decay 
at lower levels. Our quantitative model only allows for 
a single degradation rate, which leads to the intermedi- 
ate half-life of 14 min as a best-fit e value. However, the 
changing degradation rate could be rationalized under 
the assumption that the cadBA mRNA has a relatively 
weak binding affinity to the degrading enzyme, such that 
active degradation only contributes significantly at high 
mRNA levels. 

Given the compatibility of model and data, we next 
used the model to infer quantitative characteristics of 
the module that are not directly assayed. For instance, 
the activity of the central regulator CadC as a function 
of the external signals lysine, pH, and cadaverine is a 
biochemical characteristic that is pivotal to the function 
of the module, but difficult to measure directly. In our 
model, this quantitative characteristic is represented by 
the signal integration function, Eq. ([T]). On the other 
hand, the experimental data shown in Fig. |3] probes the 
final output of the Cad module on the protein level (spe- 
cific CadA activity), which is the integrated result of a 
nonlinear dynamical system with feedback. Fig. |6] com- 
pares this final system output (black curves, as in Fig. 3j 
but normalized to one) with the inferred activity of the 
CadC regulator (blue curves). The latter represent the 
three sigmoidal functions that make up the signal in- 
tegration function ([T]), with the parameters determined 
from the global fitting procedure described above. We 
observe from Fig. [6] that the final system output behaves 
qualitatively similar to the inferred biochemical activity 
of CadC. However, in each case, the apparent activation 
threshold for the system response (point of half-maximal 



CadA activity, black curve) is shifted with respect to the 
inferred biochemical activation threshold (point of half- 
maximal CadC activity, blue curve). These shifts are due 
to the fact that the total CadA activity does not only de- 
pend on the characteristics of the regulator, but also on 
the biochemical properties and timescales of the negative 
feedback loop. In principle, the feedback can even lead 
to non-monotonic behavior in the dose-response, despite 
the underlying monotonic dependence of the signal in- 
tegration function on the levels of the external signals^. 
The inferred in vivo biochemical activation thresholds, 
Ki = 3.6 mM and Kc = 235 /iM, and the values of the 
Hill coefficients describing the sensitivities to the signals, 
can be read off directly from the blue curves, as indicated 
in Fig. [6) 

It is useful to compare these biochemical activation 
thresholds to the actual concentrations encountered in 
our induction kinetics experiments, in particular at the 
time of transcriptional down-regulation {t = 30 min). 
From the plots in Fig.[4](c) and (d) we see that the lysine 
concentration at this time point is still about a factor of 3 
higher than the activation threshold, whereas the cadav- 
erine concentration of 300 jaM exceeds the deactivation 
threshold. The individual regulatory contributions from 
lysine and cadaverine to the CadC activity, i.e. g{l{t)) 
and h{c(t)) are plotted in Fig. [4] (e). The lysine curve 
(dashed line) displays only a very weak impact on the 
CadC activity, whereas the increase of cadaverine is the 
primary effect causing the down-regulation of the CadC 
activity (solid line). Hence, the analysis with our quan- 
titative model strongly suggests that the negative feed- 
back via external cadaverine can quantitatively explain 



^ The model predicts indeed a weak non-monotonic effect in the 
lysine dependence, however this is not a robust prediction. 
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FIG. 5: Correlations between the goodness of fit and the 
estimated parameters. 

The points correspond to local optima in the parameter space, 
for which the difference between the quantitative model and 
the experimental data in Figs. [2] |3] and Table [l] is minimized, 
see Materials and Methods for all details. As the fit quality 
increases (lower x^) most parameters are confined to narrow 
intervals, indicating that their values are well constrained by 
the experimental data. However, some parameters display 
significant variation even for the lowest values and from 
parameter-parameter correlation analysis in Fig. SI one finds, 
that only combinations of those are well confined by our data. 



the timing of the transient response in the wild-type Cad 
system and that the decreasing lysine stimulus is not in- 
volved in this behavior. 



Prediction and experimental analysis under altered 
conditions 



So far we have analyzed the Cad module only in the 
wild-type strain and only with a single induction proto- 
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FIG. 6: Comparison of the extracted signal integration func- 
tions {blue curves) with the normalized dose-response curves 
(black curves). 

The signal integration functions describe the dependency of 
the CadC activity on the pH level (a), on lysine (b) and on 
cadaverine (c). They correspond to the functions /, g and h 
in Eqs. ^ and pf and are plotted for the best fit parameters 
listed in Table The apparent activation thresholds of the 
dose-response curves (black curves) are shifted with respect 
to the inferred biochemical thresholds of CadC (blue curves), 
since the dose-response curves also depend on the biochemi- 
cal properties of the nonlinear feedback in the Cad module. 
From the extracted signal integration functions we can also 
read off the sensitivity of CadC on its input signals, since their 
maximal slopes are determined by the Hill coefficients ni and 



col. To obtain a more complete picture of its quantita- 
tive behavior, we constructed a mutant strain MG1655- 
lysP211 with a truncated and inactive form of the lysine 
permease LysP, see Materials and Methods for details. 
Qualitatively, we expected that this mutation would com- 
pletely abolish the lysine requirement for the activation 
of CadC, since the truncated form of LysP would be un- 
able to repress CadC [21 (an early study also indicated 
a de-repressed activation of the Cad system by a spon- 
taneous mutation in lysP [9^). Within our quantitative 
model, the lysP211 mutation was mimicked by setting 
the lysine-dependent activity function in Eq.(T]) equal to 
its maximal value, i.e., g(l) = 1, thereby rendering it in- 
dependent of the inducing lysine concentration. Apart 
from this "m silico mutation" , we left the model and the 
parameter values unchanged. 

We verified that the LysP-deficient mutant was indeed 
inducible by a shift from pH 7.6 to 5.8 alone, and did not 
require lysine for its induction (data not shown). Then, 
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we performed the same induction kinetics experiments 
as with the wild- type strain, see the results in Fig. [4] for 
the lysP211 mutant (red triangles) and the quantitative 
model (red curves). We clearly see that the expression 
in the mutant remains transient, supporting again the 
conclusion that a signal different from lysine mediates 
the transcriptional down-regulation. The shape of the 
response is altered, however, with a stronger initial ex- 
pression and a peak in the mRNA level that has a larger 
amplitude and reaches its maximum earlier than for the 
wild-type strain. Exactly these features are expected also 
on the basis of the quantitative model: The stronger ini- 
tial expression is due to the full relief of the LysP repres- 
sion, and the negative feedback via cadaverine sets in 
at an earlier time since the cadaverine threshold is more 
rapidly reached, see Figs. El (h) and (j). Also, the Cad A 
activity, shown in Fig.|4](g7, is expected to reach a higher 
steady state plateau, as observed in the experiment. In- 
terestingly, at the end of the experiment, the cadaverine 
level in Fig. [4] (h) reaches almost the 10 mM level of ini- 
tially added lysine, which is in line with our flux balance 
assumption in the model (see previous section). 

Next, we considered altered environmental conditions 
for the induction of the Cad module. Neely et al. [12^ 
had already shown that induction with initially added 
cadaverine causes a significantly reduced long-term ac- 
tivity of the cad operon, cf. Table |lj Yet, it is not clear 
whether this diminished long-term activity is caused by 
a cadBA expression pulse of similar strength but reduced 
duration, by a pulse of reduced strength with similar du- 
ration, or by a combination of both. To resolve this ques- 
tion, we again performed kinetic induction experiments 
with the wild-type strain, under identical conditions as 
in Fig. [2) but additionally with 80 /iM or 320 /iM cadav- 
erine supplied at the time of induction. In the latter case, 
the initially supplied cadaverine already slightly exceeded 
the inferred inactivation threshold of Kc = 235 /iM, such 
that a strong effect on the response could be expected. 
The resulting data is shown in Fig. [t] (a) (squares and 
triangles), together with the original data (no added ca- 
daverine, circles) for comparison. Fig. [t] (b) shows the 
same data, but with all curves normalized to peak height 
1, in order to emphasize the shape of the response. We 
observe that the primary effect of the addition of initial 
cadaverine is to reduce the strength of the response. This 
is also predicted by the quantitative model (solid lines), 
rather accurately for the 80 /iM cadaverine data set, while 
the reduction for 320 /iM cadaverine is predicted to be 
stronger than observed experimentally. 

It should be noted, that transcription was detected 
by Northern blot analysis, so that the shape of the re- 
sponse could be more accurately determined than the 
absolute amplitude. It is then interesting to observe that 
the quantitative model predicts a widening of the expres- 
sion peak and a shift of the maximum to a later time for 
the highest cadaverine concentration (320 /iM, red curve). 
Again, this behavior results from the nonlinear feedback 
in the Cad module, i.e. the addition of initial cadaverine 




Time after induction [min] 

FIG. 7: Experimental test of the kinetic model prediction. 
The lines show the parameter-free model prediction for induc- 
tion of the Cad system by a shift to pH 5.8 and 10 mM lysine, 
together with the cadaverine concentration indicated in the 
legend. The experimental data in (a) (symbols) was recorded 
and scaled as described in Materials and Methods. In (b) all 
data were normalized to their maximal values. 



decreases the initial rate of mRNA production such that 
the CadA level increases more slowly, and the negative 
feedback also sets in more slowly. This predicted change 
in the shape of the expression peak for 320 /iM agrees 
remarkably well with the experimental observation, see 
Fig. [t] (b) . This finding provides strong evidence that the 
quantitative characteristics of the Cad module are well 
described and understood with the help of our mathe- 
matical model. 



Discussion 

Conditional stress response with feedback inhibition 

In this work we analyzed the kinetics of a conditional 
pH stress response system, the lysine-decarboxylase sys- 
tem of E. coli, which exhibits only transient induction, 
even when the pH stress persists. Our results strongly 
suggest that the additional stimulus for the conditional 
response, i.e. a lysine-rich environment, is also not re- 
sponsible for the transient behavior. Rather, our kinetic 
and dose-response experiments in combination with our 
quantitative model clearly indicate that a negative feed- 
back via the product of the decarboxylation reaction, ca- 
daverine, leads to the down-regulation of the response. 

Cadaverine has previously been linked to the transient 
behavior, in pivotal work on this stress response sys- 
tem [12j [13]. This link was based on the observations 
that external addition of cadaverine significantly reduces 
the long-term Cad activity and that a CadA~ mutant 
displayed persistent cadBA expression. However, if the 
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transcriptional shut-off was mediated by a decrease of the 
lysine stimulus, the persistent cadBA expression could 
alternatively be explained by the lack of lysine consump- 
tion in the CadA~ mutant. And while a reduction in 
the steady-state activity implies that the Cad system is 
generally repressed by cadaverine, it is not clear that 
the timing of its down-regulation is set via this nega- 
tive regulatory interaction. None of the previous studies 
directly measured the system-induced dynamics of the 
cadaverine concentration or studied the kinetics of the 
Cad system with externally added cadaverine. By per- 
forming these quantitative experiments with a high time 
resolution and by interpreting them with a quantitative 
theoretical model, we obtained evidence for a causal rela- 
tion between the time of transcriptional down-regulation 
and the increase of the external cadaverine concentra- 
tion above its deactivation threshold. For instance, the 
decreasing amplitude in the dynamical response after 
adding external cadaverine, as well as a more subtle de- 
layed down-regulation, both predicted by the quantita- 
tive model, were strikingly confirmed by our kinetic mea- 
surements. 



Signal transduction mechanism 

Our systems-level study of the Cad module also per- 
mits some conclusions about the involved molecular in- 
teractions and the signal transduction mechanism. Us- 
ing our quantitative model, we were able to estimate the 
relevant in vivo activation and deactivation thresholds 
of the regulatory protein CadC. For the inactivation of 
CadC by external cadaverine we found a threshold of 
Kc = 235 /iM. This value is surprisingly close to the in 
vitro binding constant of 96 /iM for the interaction of ca- 
daverine with the periplasmic domain of CadC [21,. In 
contrast, CadC has almost no affinity for lysine. There 
is recent evidence that CadC is inhibited at low lysine 
concentrations via a transmembrane domain interaction 
with the lysine permease LysP, whereas the interaction is 
released at high lysine levels [21] . In the present work we 
determined the effective in vivo lysine activation thresh- 
old to be i^^ = 3.6 mM. This result is somewhat surpris- 
ing, since the Michaelis constant Km for lysine trans- 
port by LysP is much lower at ~ 10 /iM [30]. However, 
the LysP/CadC interaction and the Km of LysP do not 
necessarily need to have a direct correspondence. 

An interesting open question concerns the signal trans- 
duction mechanism of CadC. Two alternative models 
have been proposed: (i) A reversible conformational tran- 
sition of CadC activates its cytoplasmic N-terminal do- 
main and allows it to bind to the promoter while remain- 
ing integrated in the membrane, (ii) A shift to acidic 
pH induces cleavage of the cytoplasmic domain, allowing 
it to diffuse freely to the promoter. The existence of the 
negative feedback by external cadaverine puts causal con- 
straints on these microscopic mechanisms. If the cleavage 
mechanism were realized, it is not clear how external ca- 



daverine could down-regulate cadBA expression after it 
has been induced, since the freely-diffusing cytoplasmic 
domains would no longer be able to recognize this sig- 
nal. A high turnover of cleaved CadC could solve this 
problem, by rapidly eliminating unresponsive activators. 
However, the increased degradation of CadC under in- 
ducing conditions would have to be balanced by an ele- 
vated cadC expression. This, however, is in disagreement 
with previous observations, where it was found that cadC 
expression is constitutive [161 [17]. Hence, the existence 
of the negative feedback by external cadaverine leads us 
to favor the reversible model. 



Top-down system analysis of a functional module 

The Cad system of E. coli is a functional module with 
few closely connected molecular components. How such 
modules integrate, process, and respond to external sig- 
nals is a central question, but generally also a difficult 
one. The approach taken in the current study is akin to a 
"top-down" system analysis, where input signals are con- 
trolled and the output (s), as well as key internal system 
variables, are measured. This is in contrast to a biochem- 
ical "bottom-up" approach, where each component would 
first be characterized separately and then their pairwise 
interactions, gradually moving upwards in complexity. In 
our study, we were able to bring these two complemen- 
tary approaches into a first contact for the Cad module, 
with the help of our quantitative model which provided 
the means to estimate relevant in vivo values for biochem- 
ical interaction parameters and the quantitative form of 
the signal integration function displayed in Fig. [6] With- 
out the quantitative model, we would have been unable 
to extract this "hidden information" from the experimen- 
tal data. Direct in vivo measurements of signal integra- 
tion functions have been performed for the regulatory 
circuit controling chemotaxis in E. coli, using sophisti- 
cated single-molecule techniques. Since these powerful 
techniques are not easily transferred to the large class 
of functional modules of interest, the indirect approach 
taken in the present work may often be a welcome alter- 
native. 



Conclusion and Outlook 

Our quantitative analysis of the Cad system provides 
a first step towards a system-level understanding of the 
complex acid stress response network in microbes. Anal- 
ogous quantitative studies of the two other major amino- 
acid decarboxylase systems, glutamate and arginine de- 
carboxylase [8 , could reveal important insights in how 
these modules are orchestrated in the complex environ- 
ment of its host. The presence of multiple amino acids 
in the natural environment of E. coli suggests that these 
conditional stress-response systems are often induced in 
parallel. It will be interesting to study how these systems 
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are coordinated to provide an effective and robust way of 
pH homeostasis and acid tolerance response. 



Materials and Methods 

Bacterial strains and growth conditions. E. coli 
MG1655 [31 was used as wild- type strain. The lysP211 
mutant was obtained via undirected mutagenesis P [32]. 
For this purpose cells of E. coli MG1655 were grown 
on minimal agar plates with 0.2% (w/v) glucose as sole 
carbon source [33] containing 100 /ig/ml thiosine (S- 
aminoethyl cysteine) , a toxic lysine analog which leads to 
spontaneous mutations in lysP. One mutant, designated 
MG1655-/^5P^77, had a nucleotide exchange at position 
211 in lysP resulting in a stop codon and hence a trun- 
cated and inactive form of LysP (70 amino acids). E. coli 
strains MG1655 and M.Gl^bb-lysP211 were grown aero- 
bically in shaking flasks under non-inducing conditions at 
pH 7.6 in 5 L phosphate buffered minimal medium [34] 
containing 0.4% (w/v) glucose as sole carbon source at 
37°C to an ODeoo = 0.5 (non-inducing conditions). Sub- 
sequently, cells were collected by centrifugation (10 min, 
4000 g at 37°C), and transferred into fresh prewarmed 
minimal medium, pH 5.8 containing 10 mM L- lysine (L- 
lysine- hydrochloride. Roth) and 0.4% (w/v) glucose in a 
5 L fermenter (Biostat B, Satorius BBI Systems GmbH) 
(inducing conditions). Cultivation of cells was contin- 
ued anaerobically at 37°C. At the indicated times, sam- 
ples were taken, centrifuged at 4000 g (4°C) for 5 min, 
and cell pellets as well as supernatants were separately 
stored at -80°C until further use. The number of colony 
forming units was determined after incubation of 100 /iL 
of various dilutions on LB agar plates overnight at 37°C 

Lysine Decarboxylase Assay. Speciflc activity of the 
lysine decarboxylase CadA was measured by resuspend- 
ing cells corresponding to 10 mL culture in 1 mL Ldc 
buffer (100 mM Na-acetate, pH 6.0, 1 mM EDTA, 0.1 
mM pyridoxal phosphate, 10 mM /3-mercaptoethanol, 
10% (w/v) glycerol). Lysozyme to 0.1 mg/mL was added, 
and the mixture was incubated on ice for 30 min. The cell 
suspension was sonifled (10-20% amplitude; 0.5 sec out- 
put puis; digital Branson Sonifler II 250), and the lysate 
was centrifuged at 15000 g at 4°C for 15 min. Activity of 
lysine decarboxylase in cell-free extracts was measured 
as described [36 using 5 jig protein per assay. Speciflc 
activity is deflned as lU/mg= 1 /imol cadaverine/(min 
X mg protein). 

Measurement of extracellular cadaverine. The extra- 
cellular cadaverine concentration was determined accord- 
ing to the spectrophotometric method described by Phan 
et al. [37 . Briefly, 10 /iL of culture supernatant was di- 
luted 5 fold with HsOdest, then 120 /iL NasCOs (1 M) 
and 120 /iL TNBS (2,4,6-trinitrobenzene-sulfonic acid; 
10 mM) were added, and the mixture was incubated 
for 4 min at 40°C. After extraction with 1 mL toluene, 
the absorption of the organic phase (containing N,N- 



bistrinitrophenyl-cadaverine) at 340 nm was measured. 
The cadaverine concentration was calculated based on a 
standard curve using cadaverine-dihydrochloride (Sigma) 
between and 500 nmol. 

Preparation of RNA. Total RNA was isolated ac- 
cording to the method of Aiba et al. [38]. Briefly, 
cells were resuspended in cold 20 mM Tris/HCl, pH 8.0, 
and subsequently lysed by addition of 20 mM sodium 
acetate, pH 5.5, 0.5% (w/v) SDS, 1 mM EDTA, pH 
8.0. Then, RNA was extracted with prewarmed (60°C) 
acid phenol, and the mixture was centrifuged at 12000 
g. After an additional extraction of RNA using phe- 
nol/chloroform/isoamylalcohol (25:24:1), RNA was pre- 
cipitated with 100% ethanol at -20°C overnight. The 
precipitate was washed with 70% (v/v) ethanol, and the 
dry RNA pellet was dissolved in 35 /iL H2 0dest- RNA 
concentration was determined by measuring the absorp- 
tion at 260 nm. All solutions were prepared with 0.1% 
(v/v) DEPC (diethylpyrocarbonate). 

Northern blot analysis. Northern blot analysis was 
performed following the protocol described earlier [39j. 
Briefly, 20 jig RNA was separated by electrophoresis 
in 1.2% (w/v) agarose- 1.1% (v/v) formaldehyde gels in 
MOPS (morpholinepropanesulfonic acid) buffer. RNA 
was transferred to a Hybond-Nylon membrane (GE 
Healthcare) by capillar blotting. Hybridization was per- 
formed following a standard protocol [35 using a [a^^- 
P]dCTP labelled PGR fragment of the first 400 bp of 
cadBA. Radioactive labelling was quantified with a Phos- 
phoimager. As control, expression of rpoD^ a house- 
keeping gene of E. coli was analyzed. Signal intensity 
of cadBA mRNA was normalized to the signal intensity 
of rpoD mRNA. If not indicated otherwise, the data are 
given as fold-change of cadBA transcription relative to 
the pre-induction value. Additionally, in Figs. |4] and [T] 
the absolute magnitude of the mRNA fold-change was 
rescaled, such that the integral over the expression curve 
was proportional to the long-term GadA activity. 

Model details. From thermodynamic models of tran- 
scriptional regulation reviewed in [26[ [27] , the effective 
transcription rate, ^ as a function of two activators A 
and B with independent binding sites (binding constants 
Ka and i^T^) is given by 



, + {A/KA)fA 

1 + [A/Ka) 



■ + {B/KB)fB 

1 + {B/Kb) 



(8) 



cf. ref. [26] Table 1, Case 10. We make the simplify- 
ing assumption that both binding sites are identical and 
that each site can only be bound by a dimer of GadG 
(C'2) JB]. Setting A = B = €2^ exploiting mass action 
K = IC2 and introducing the effective binding con- 



stant Kc = ^/KKc2 leads to the first term in Eq. (|4|. 

The effective Michaelis-Menten form of the lysme 
turnover rate in Eq. ([6| was derived as follows. The trans- 
port of lysine and cadaverine via GadB was modeled in 
analogy to the homologous arginine-ornithine antiporter 
ArcD in L. lactis [40], see Fig.jsJ For low external cadav- 
erine and internal lysine concentrations a general form 
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the inwardly directed flux of lysine is [40] 



^eff 



= ktB- 



K 



Bl 



(9) 



where B is the number of CadB molecules per cell and 
[/] and [c] are the concentrations of lysine and cadaverine 
on the outer surface of the membrane, respectively. In- 
ternal solute concentrations are marked by a prime C). 
The parameters ^ Kbi and Kbc' are deflned in the 
kinetic scheme of Fig. [sj Here it is assumed that (i) the 
conformational transition in CadB, which mediates the 
transport, does not occur without bound lysine or cadav- 
erine, (ii) the antiporter and its substrates are in binding 
equilibrium at each surface of the membrane, and (iii) 
the membrane translocation reaction of the carrier is slow 
and rate-limiting. The internal cadaverine concentration 
[c'] is determined by an interplay of the (reversible) de- 
carboxylation through CadA pT and by the export via 
CadB. In steady state we flnd 



A A 



(10) 



where A is the number of CadA molecules per cell. If 
the equilibration between lysine and cadaverine through 
CadA is fast compared to the transport through CadB, 
the second term is negligible and the internal cadaverine 
level is solely determined by the internal lysine level [V]. 
In steady state it turns out that [V] is not affected by 
the Cad module, since the 1:1 stoichiometry of the an- 
tiporter assures that lysine decarboxylation is balanced 
by lysine import. If we further take advantage of the 
fact that CadB and CadA are transcribed polycystroni- 
cally and that there seems to be no post-transcriptional 
regulation [11 , we can set CadB proportional to CadA, 
i.e., B = a A. Taken together, in the limit of rapid 
CadA kinetics and low external cadaverine and internal 
lysine concentrations Eq. ([9| reduces to the simple ef- 
fective Michaelis-Menten form in Eq. (|6| with CadA as 
the enzyme, Vmax = (o^ x /c^"^)/?7 as the effective maximal 
turnover rate and K = Kbi/t] the effective Michaelis 

constant, where 77 = 1 + 



kc [U] 



Parameter estimation. The parameters of our quan- 
titative model were estimated by using a trust-region 
reflective Newton method (MATLAB, The MathWorks, 
Inc.) to minimize the total x^? deflned by 

X'(^l = xlM + xIh0) + Xls0) + Xlad0) , (11) 

with respect to the parameter vector = (^1, . . . •,Om)i 
where M = 14 is the total number of model parame- 
ters, cf. Table [Til The contribution of the kinetic data is 
calculated from 
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FIG. 8: Kinetic scheme of lysine and cadaverine transport 
and turnover. Antiport of lysine (/) and cadaverine (c) by 
CadB is modeled in analogy to the homologous arginine- 
ornithine antiporter ArcD in L. lactis by a single site Ping 
Pong Bi-Bi mechanism [40^. Interconversion of lysine and 
cadaverine by CadA is modeled by a reversible first order 
reaction [JJ. The prime (') indicates internal quantities. 
The rates {k^ , , , k~ , k^, k^) and equilibrium constants 
[Kbu Kbv ^ Kbc^ Kbc') are indicated next to the reaction 
steps. 



where yij^ y2j and ysj are the experimental data of the 
cadBA mRNA, the CadA activity and cadaverine concen- 
tration at time tj, respectively. Similarly, the yij{0) de- 
note the corresponding values of the quantitative model 
for a given parameter set 6>, and the ciij are the standard 
errors of each measurement (estimates from our experi- 
ments: (Jij = 5, (J2j = 0.1 U/(mg protein), and osj = 0.5 
mM Vj). The Ni are the number of datapoints of a given 
timeseries, Atij = — ti(j_i))/2 is the time be- 

tween subsequent data points and Ati is the mean time 
between the data points in dataset i. The contribution 
of the dose-response curves are similarly deflned and are 
exemplarily shown for the pH-dependent response 
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pH 



'^f Ar{pHi)-Ar{pHi,( 



i=l 



(13) 



where A^{pHi) and A^{pHi.,0) are the experimental and 
theoretical CadA activities after time r at a given pH 
\e\e\ pHi {i = l,...A/p//), respectively. The cr^ denote 
the standard errors of the measurement and NpH is the 
total number of datapoints. The other contributions xfys 
and xtad deflned equivalent ly to Eq. (13). However, 
they differ in the time at which the CadAactivity was 
determined experimentally, i.e., TpH = 1.5 h, riys = 8h, 
and Tcad = 3 h. 

To account for the presence of local optima, and to 
quantify the uncertainty in the estimated parameters, we 
performed 1000 independent flts with randomly chosen 
initial parameter sets (within their physiological ranges) . 
In Fig. [5] the flnal values are plotted against the flnal 
parameters. We followed Ref. [42] to compute the asym- 
metric errors cr+ and cr_ with respect to the optimal pa- 
rameter values O'^P^ listed in Table [nj The squared errors 
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for parameter 0]^ were calculated using the equations 



i:0k,i>6 



, and 



i-ek.i<6T 



E 



-X?/2 



(14) 



where 6]^^i is the value of parameter 6^ in the ith fit, ^^^^ 
is the value of Ok in the fit with the lowest value of x^, 



and Xi is the value of x^ for the ith fit. In using the like- 
lihood function we assume that the errors in the 
measurements are independent and normally distributed 
with widths equal to the standard error of the mean. 
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